Quantum diffusion in biased washboard potentials: strong friction limit 
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Difltusive transport properties of a quantum Brownian particle moving in a tilted spatially pe- 
riodic potential and strongly interacting with a thermostat are explored. Apart from the average 
stationary velocity, we foremost investigate the diffusive behavior by evaluating the effective diffu- 
sion coefficient together with the corresponding Peclet number. Corrections due to quantum effects, 
such as quantum tunneling and quantum fluctuations, are shown to substantially enhance the ef- 
fectiveness of diffusive transport if only the thermostat temperature resides within an appropriate 
interval of intermediate values. 
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INTRODUCTION 

Brownian motion in periodic structures can describe 
diverse processes in many different branches of science. 
Within a physical context, among other phenomena, it 
models the dynamics of the phase difference across a 
Josephsonj unction 0, Q, rotating dipoles in external 
fields y, y|, superionic conductors charge density 
waves |3 , particle separation by electrophoresis jTj , trans- 
port on crystalline surfaces [a 0, biophys ical processes 
such as intracellular transport |lCl IllL Il2l Il3| . Yet an- 
other important area constitutes the noise-assisted trans- 
port of Brownian particles 0, 0| , as it occurs for Brow- 
nian motors possessing ample applications in physics and 
chemistry 10]. 

In this paper, we study the one-dimensional over- 
damped motion of a quantum Brownian particle sub- 
jected to a tilted potential U{x), 



U{x) ^V{x) - Fx, V{x) ^V{x + L), 



(1) 



where V{x) denotes a periodic potential of period L and 
F is an external static force. 

The basic quantities characterizing this motion are sta- 
tistical moments of position and velocity of the Brownian 
particle. At least the first two moments i.e. the average 
position and average velocity and their respective disper- 
sions are most substantial. In particular, the stationary 
average velocity can be defined by the relation 



(.) = hm Ml, 

t^oo t 



(2) 



where x{t) is the position of the Brownian particle at 
time t and (...) is the average over all realizations of the 
thermal noise and the initial conditions. The dispersion 
of the position can be characterized by the diffusion co- 
efficient defined as 0| 



In the classical case, for strong friction, when only ther- 
mal equilibrium ffuctuations affect the particle, the sta- 
tionary average velocity and the diffusion coefficient can 
be expressed by exact closed formulas 0, For 
nonequilibrium dri ving the y ca n be calculated in specific 
cases only, see e.g. |l9H2dl2lj . 

Depending on the form of the potential, the magnitude 
of the tilt and thermostat temperature, two interesting 
phenomena can be observed: A giant enhancement of the 
diffusion constant at a critical tilt llVt lisll and a low ran- 
domness window at sub-critical tilts |22l |. The influence 
of the shape of the potential E3 . and of a position 
dependent friction coefficient '24] on the transport have 
also been investigated. 
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FIG. 1: Schematic picture of the Brownian particle in the 
tilted periodic potential U{x) = V{x) — Fx, defined in Q 
and 1171 1 below, plotted for the following set of parameters: 
A = 10, e = 10, _F = 0.2, L denotes the period the potential 
V{x) and Lr stands for the length of free slide between the 
barriers. 



However, in many cases, like for the Josephson junc- 
tion at intermediate temperatures, the classical the ory 
is insufficient; i.e., the leading quantum corrections |25j 
should be considered. It was shown in Ref. ^2^] that 
in the stron g fri ction limit, effects of quantum Brownian 
ffuctuations [28| are restricted not only to low tempera- 
tures; therefore, these should be incorporated for higher 
temperatures as well. This is so because quantum ffuctu- 
ations, even if reduced for one variable, are enlarged for 
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the conjugate variable. Quantum corrections can modify 
the dynamics quantitatively and sometimes even quali- 
tatively. Physically relevant exam ples illustrating these 
features are presented in Refs. |27L I2M l2g| . 

For the average current in of the Brownian particle 
dynamics such quantum corrections have been studied 
repeatedl y in the previous literature within different ap- 
proaches |30j. In particular, within the quantum Smolu- 
chowski equation these quantum corrections have been 
studied recently by Ankerhold in Ref. for an over- 
damped Josephson junction. 

In distinct contrast, with this work we mainly focus 
on the role of the quantum corrections to the diffusion 
of the mean squared displacement of the coordinate (or 
Josephson- phase, respectively) as defined with 

The quantum diffusive dynamics in the regime of 
strong interaction with a thermostat can be described by 
the so called Quantum Smoluchowski Equation (QSE) 
incorporating quantum fluctuations above the 
crossover temperature [ij, . It corresponds to a clas- 
sical Smoluchowski equation, in which the potential U{x) 
and the diffusion coefficient D are modified due to quan- 
tum effects. In other words, the quantum non-Markovian 
diffusion process of a particle position is approximated 
by a classical Markov process describing a motion in an 
effective potential and with an effective, state-dependent 
diffusion coefficient. This leads to a quite comfortable sit- 
uation because methods of analyzing Smoluchowski equa- 
tions are well elaborated and can directly be applied and 
implemented for a wide class of systems. In this paper we 
demonstrate how the quantum fluctuations influence the 
diffusion behavior of a Brownian particle stochastically 
moving in washboard-like potentials. 



CLASSICAL BROWNIAN PARTICLE IN TILTED 
PERIODIC POTENTIALS 

The overdamped motion of a classical Brownian parti- 
cle is described by the Langevin equation 



■qx = -U'{x) + y%kT £,{t), 



(4) 



where rj denotes the viscous friction coefficient and k is 
the Boltzmann constant. The dot and the prime indicate 
differentiation with respect to time t and to position x, 
respectively. The zero-mean and 5-correlated Gaussian 
white noise ^(t) models the influence of a thermostat of 
temperature T on the system. 

The calculation of the average velocity (v) and effec- 
tive diffusion D^f / can be accomplished by mapping the 
washboard potential on a corresponding jump process. 
This construction procedure has been elucidated in [2^ . 
As a result, a cumulative process with independent incre- 
ments is obtained and its asymptotic mean and variance 
are given by the first two central moments of the escape 



time density IT3.|l7l ll^. In this way, the stationary av- 
era ge y elocitv |l4lll5ll32L Is^ and the diffusion coefficient 
[iTlllSl l for the process modelled by Q can be expressed 
by closed form relations involving quadratures only; i.e., 

L 



{v) 



D. 



Ti{xo XQ + L)' 
2 



(5) 



T?ixo^xo + L) 



Tl{xo^xo + L) 
where xq is an arbitrary, initial value and 
Tnixo^b) = {r{xo^b)) 



(6) 



denotes the nth statistical moment of the first passage 
time t(xo — > b) at which the Brownian particle arrives 
at the point b while starting out from the position xq. 
For the case b > xn, these moments are given by the 
recurrence relation 1361 . 



Tn{xo b) — nfirj I dx exp[/3[/(a;)] 



dy exp[-/3[/(y)]T„_i(2/^6) 



(7) 



for n = 1,2,3..., where Ta{y b) = 1, (3 = 1/kT and 
the product (3?] — D^^^ is the reciprocal of the Einstein 
diffusion coefficient Dq. The expressions © and {Tj) are 
rather complicated. However, they can be simplified as 
shown in [l7llli|. 



OVERDAMPED QUANTUM BROWNIAN 
MOTION 



To start with the investigation of quantum corrections 
to diffusion we consider a quantum Brownian particle 
moving in the tilted potential U{x). The evolution of its 
position can be described by the respective probability 
density function P{x,t) = {x\p(t)\x), which is the di- 
agonal part of the statistical operator p{t). Within the 
strong friction limit (the quantum Smoluchowski regime) , 
the dynamics of such a particle is described by the Quan- 
tum Smoluchowski Equation (QSE) that takes into ac- 
count leading quantum corrections. It has the structure 
of a classical Smoluchowski equation with modified drift 
and modified diffusion terms 12611271 l29l 



The effective potential reads 

Ueff{x)^U{x) + {l/2)XU"ix). (9) 

The effective diffusion coefficient D{x), being constant in 
the classical case, i.e., D{x) = D = ksT = /?~^, becomes 
position-dependent, assuming the unique form [29ll34|. 



d 



d 



D{x) = {l3[l-\f}U"{x)] 



(10) 
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This diffusion is required to remain non-negative, i.e., 
within its regime of vahdity [29I Ib^ . the inequahty 
X/3U"{x) = \pV"{x) < 1 must be satisfied for all po- 
sitions X. For smooth periodic functions V{x) and suffi- 
ciently small A/3 this inequality holds for arbitrary x. 

The prominent parameter A characterizes quantum 
fluctuations in position space; it explicitly reads [1^113 j 



A = {h/Trr/) \ii{h(3r]/2TrM) 



(11) 



It depends nonlinearly on the Planck constant h and on 
the mass M of the Brownian particle, whereas, in the 
classical case, the overdamped dynamics does neither de- 
pend on h nor on the mass M (note that we use the 
friction constant rj which has the unit [kg/s] as in the 
classical Stokes case). Note also that this quantum cor- 
rection approaches zero with the friction 77 growing to- 
wards infinity. 

The Langevin equation corresponding to the 
Smoluchowski equation (jSj) becomes within the Ito- 
intcrprctation 35], 



i^x^-U'Ax) + ^2riD{x) m ■ 



(12) 



The average stationary velocity (v) and the diffusion 
classical coefficient D^f / can be calculated as in the case 
described by the Langevin equation using the rela- 
tions © and the known formula for statistical moments 
of the first passage time. In comparison with {Tj), the 
statistical moments are thereby modified into the form 

M 

Tn{xo ^ b) = nrj dx exp[(j){x)] 

J Xo 

xf dyD-\y)eM-Hy)]Tn-i{y^b), (13) 



where 



0(x) = / dz. 



(14) 



Insertion of the expressions for the effective potential © 
and the effective diffusion function IjlOl) yields 

Tn{xo ^ b) — n(3ri / dx exp[/3'0(a;)] 



X / dy exp[-my)] [1 - XpU"iy)] T^-i{y ^ 6),(15) 

J —QO 

where the thermodynamic potential ip{x) becomes 

V'(x) = Uix) + {l/2)XU"ix)+ (16) 
- {l/2)Xp[U'{x)f-{l/4)X^(3[U"{x)f. 



classical form (|7I) in a two-fold way: First, the physical 
potential U (x) is replaced by the thermodynamic poten- 
tial 'ijj{x). This thermodynamic potential depends on the 
temperature /? of the system and on the coupling con- 
stant of the Brownian particle with its surroundings via 
the damping constant 77, which in turn enters into the 
parameter A. Second, the function in the inner integral 
(over the variable y) on the right hand side of eq. H15|l 
is modified by the factor [1 — XpU"{y)], which depends 
on the curvature (i.e. on U"{y) = V"{y)) of the physical 
potential U{y). 

We will examine the quantum Brownian particle mov- 
ing in a tilted washboard potential like that presented in 
Fig. n We scale the force and the diffusion coefficient in 
such a way that for the rescaled equation corresponding 
to (|12|l the friction coefficient equals rj ~ 1. We choose 
a specific form of the periodic part V{x) of the potential 
U(x), namely 



V{x) = A exp[£(cos(x) — l)]/e 



(17) 



The advantage of this choice is that by an appropriate 
manipulation of e and A, the barrier height and the dis- 
tance between neighboring barriers can be varied inde- 
pendently. As a consequence, one can change two time- 
scales independently: a first one is related to the deter- 
ministic sliding motion between neighboring barriers and 
the other is related to the inverse of the activation rate 
over a barrier. 



QUANTUM DIFFUSION IN TILTED 
WASHBOARD POTENTIALS 



In order to understand the influence of the shape of 
the washboard potential on the particle dynamics, both 
in the classical and quantum regimes, is it desirable to 
identify characteristic time-scales. The first one is given 
by the time — Lr/F the particle needs to slide down 
the distance Lr (see in Fignj with a constant velocity 
V = F (remember that the friction coefficient 77 = 1 and 
the force is rescaled). This timescale is relevant if the po- 
tential has an almost constant slope between neighboring 
comparatively narrow barriers, as in the case considered 
here. The second time-scale Te is determined by the es- 
cape time over the barrier. The potential V{x) has been 
chosen in the above described way, so that these time- 
scales can be 'tuned' independently: by the force F 
and the parameter e, and by the barrier height using 
both parameters e and A of the potential (|17|l . 

The transport of particles is optimal if a large mean 
velocity goes along with small diffusion. This can be 
quantified by the dimensionless Peclet number j37j : 



We observe that quantum corrections modify the sta- 
tistical moments as given by eq. (|15|) compared to the 



(18) 
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FIG. 2: Average velocity (v) (a), effective diffusion -De// (b) 
and Peclet number Pe (c) are drawn versus tilt F. Parame- 
ters of the potential are: e = 100, A = 10, which result in a 
barrier height of AV = 0.1 in the absence of a tilt {F — 0), 
and a critical tilt Fc — 0.6 above which all barriers disappear. 
The rescaled inverse temperature is set to l3 — 100, corre- 
sponding to a small quantum parameter Ao — 0.00115129. 
The mean velocity and Peclet number are always larger in a 
quantum case, even if the classical effective diffusion is smaller 
(as is the case for F < 0.25). With the panel (d) we depict 
the corresponding classical, U{x), (dashed), quantum effec- 
tive, Ueff{x), (solid line) and thermodynamic, ip{x), (dotted) 
tilted potentials, respectively. The position-dependent diffu- 
sion coefficient is plotted in panel (e). 
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FIG. 3: Average velocity (a), effective diffusion (b) and Peclet 
number (c) versus the inverse temperature /3. The chosen pa- 
rameters are: bias F = 0.2, potential parameters e = 100 and 
A — 10. The two arrows in panel (c) mark the characteris- 
tic temperatures where the system starts to markedly deviate 
from the classical behavior (see text for details). The region of 
intermediate temperatures is located between the downward 
and upward pointing arrows. It bridges between the regions 
of freely sliding Brownian particle dynamics, characterized by 
the Peclet number Pe — FL[3 (see the thin, dotted line in (c)) 
and the Poisson like behavior with Pe = 2. Within this region 
the Peclet number assumes its maximal value. 



The efficiency of the diffusive transport as measured by 
the Peclet number can either be enhanced by an increase 
of the net current (i.e. the stationary mean velocity) 
and/or by a decrease of the effective diffusion, resulting 
in a maximal Peclet number Pe. The average velocity for 
the overdamped motion in a tilted washboard potential 
is limited by the free-slide speed which coincides with 
the value F in our case. The particle will approach this 
free-slide velocity when the barriers become negligible. 
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for example for a sufficiently high thermal energy kT or 
very strong tilt F. This situation, however, does not 
lead to an optimal transport performance in the sense of 
a maximal Peclet number 1221 1381. see also below. 
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FIG. 4: The influence of quantum corrections are illustrated 
as a function of bias F vs. inverse temperature /3 by the rela- 
tive difference between the quantum (Q) and the correspond- 
ing classical (C) values of the current {{v)'^ — (v)'^) / {v}^ (a), 



the effective diffusion {D^^f 



Deff)/F>^ff (b) and the Peclet 



Q 

e/. 

number (Pe*^ — Pe'^)/ Pe^' (c), respectively. The positive 
values of the depicted quantities indicate the relevant role of 
quantum effects. 



Quantum-renormalization of the barrier shape 

We studied the influence of quantum corrections on 
transport in tilted periodic systems by means of a numer- 



ical analysis of the basic expressions ® and lfTB|) . We 
found that the quantum current is always higher than 
the corresponding classical one (see Figs. |21 El and 0)) . 
This phenomenon can be explained by comparing the 
potential U{x), the effective quantum potential Ueff{x) 
and the thermodynamic potential ip{x) with each other, 
as well as by analyzing the effective diffusion function 
D{x) (which is constant in the classical case), see Fig. |21 
Clearly, the effective quantum potential Ueff{x) (solid 
line) possesses slightly lower and thinner barriers than 
U{x) (dashed hue). 

The state-dependent diffusion function D{x) possesses 
maxima and minima. The maxima, which are shifted 
away from the potential barrier locations, can be inter- 
preted as a higher effective local temperature. The mini- 
mum of D{x) is located in the neighborhood of the top of 
the barrier (near a; ~ in panel (e)). ft means that quan- 
tum fluctuations mimic an effective temperature which 
is lower at the barrier and higher in the potential wells. 
For the escape dynamics the thermodynamic potential 
ipix) is decisive: It contains the combined influences of 
the effective potential and the effective diffusion. In the 
present case, ipi^) displays both a lower and a narrower 
barrier than the bare potential U{x) of the corresponding 
classical process. It is remarkable that for all cases con- 
sidered, the Peclet number is always larger in the quan- 
tum case than in the classical case, thus providing a more 
coherent motion. This behavior is exemplified in Figs. |21 
El and H 



Role of quantum corrections for diffusion 

Depending on the relation between the thermal energy 
kT and the barrier height of the tilted potential, one 
can distinguish regimes of high and low temperatures, 
denoted in Fig. |3|by a downward and an upward arrow, 
respectively. To the left of downward arrow, i.e. for high 
temperatures, the particle barely feels the barriers and 
thus freely slides down the hill with an average velocity 
given by (v) w F. 

The effective diffusion coefficient D^f / can then be ap- 
proximated by the Einstein diffusion coefficient in 
both the classical and in the quantum cases. The cor- 
responding approximated values of the Peclet number, 
given by Pe = FLf3, are depicted as a dotted line. To 
the right of the upward arrow in Fig. Et the relaxation 
time Tr is much smaller than the time scale for barrier 
crossing. The time evolution in that case consists of a se- 
quence of independent activations. Indeed, the value of 
Peclet number approaches Pe = 2 which is characteristic 
for Poisson process. 

The most interesting region is located between the two 
temperatures, indicated by two oppositely directed ar- 
rows (one downwards and one upwards) , where we found 
the optimal transport, i.e. the maximum of the Peclet 
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Quantum Classical 




FIG. 5: The influence of the quantum corrections on the 
evolution of the coarse grained probability density function 
P{x,t) (see discussion in 22]) is illustrated for the same set 
of parameters as in Fig. |21 and for different forces F = 0.2 
and 0.3. Note that in the quantum case the velocity is always 
larger, therefore, the maximum of the density is located at a 
more distant position as compared to the classical case. The 
width of the density differs depending on the driving param- 
eters and thereby reflects the effective diffusion strength. 



number Pe. The quantum behavior significantly devi- 
ates from the classical one only within this very region. 
We observe that the transport quality, expressed in terms 
of the characterize! Pe is never suppressed by quantum 
effects (see Fig. Oland^, even though quantum correc- 
tions may increase the the effective diffusion -De//- 

In Fig. 21 we depict the relative value for the correc- 
tion of Deff and Pe, respectively, in the parameter space 
spanned by /3, F. First, it is detectable that the corrected 
values of velocity, the effective diffusion and Peclet num- 
ber may differ up to 200% from the classical values. Sec- 
ond, the sign of the relative corrections of (w) and Pe is 
positive, but in the case of the effective diffusion De// it 
might assume negative values within some range of pa- 
rameters. 

Finally, we illustrate the effect of quantum corrections 
on the transport in tilted washboard potential. In Fig. 
we present two examples of the time evolution of the den- 
sity function P{x, t). The impact of quantum corrections 
is clearly visible. The width of P{x,t) becomes larger 
at F = 0.2 and smaller at _F = 0.3 when the quantum 
corrections are acting, but the peak of the probability 
density travels in the quantum case with a larger veloc- 
ity in both situations. 

CONCLUSIONS 

The effect of the quantum contribution on Brownian 
motion, in particular on the diffusion of particles and 
the related transport performance is addressed in this 
work. The quantum current always exceeds the corre- 
sponding classical one; - quantum features like tunnel- 



ing and quantum fluctuations seemingly always assist the 
particle to overcome barriers and to pass longer distances, 
resulting in larger average stationary velocity. The dif- 
fusion coefficient, D^ff, is found to assume, generally, 
a non-monotonic function of the static force F and the 
temperature (3. In other words - optimal conditions ex- 
ist for both, directed and diffusive transport. Depending 
on the parameters of the system, quantum effects may 
either increase or decrease the effective diffusion of the 
particle. The Peclet number is found to be always larger 
for quantum systems, - see in Fig. ^ (c). The most sig- 
nificant finding is that quantum effects always improve 
the diffusive quantum transport for the class of systems 
considered in this work. 
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